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Abstract 

We describe AlterBBN, a public C program for evaluating the abundances of the ele- 
ments generated by Big-Bang nucleosynthesis (BBN). This program enables the user 
y—{ to compute the abundances of the elements in the standard model of cosmology, and 

additionally provides possibilities to alter the assumptions of the cosmological model in 
order to study their consequences on the abundances of the elements. In particular the 
y—{ baryon-to-photon ratio and the effective number of neutrinos, as well as the expansion 

rate and the entropy content of the Universe during BBN can be modified in AlterBBN. 
Such features allow the user to test the cosmological models by confronting them to 
BBN constraints. A presentation of the physics of BBN and the features of AlterBBN 
is provided here under the form of a manual. 
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1 Introduction 



Big Bang Nucleosynthesis (BBN) is an important part of the Cosmological Standard Model, 
as it gave birth to the first nuclei. Based on BBN models it is possible to calculate the abun- 
dances of the light elements in the Universe, which can then be compared to the measured 
abundances. This allows to probe the Universe properties at the time of nucleosynthesis, 
which is the most ancient period observationally accessible. In standard cosmology, only 
one free parameter, the baryon-to-photon ratio at BBN time is needed to determine the 
abundances of the elements, and this parameter can be deduced from the observation of 
the Cosmic Microwave Background pQ . To calculate the abundances of the elements in the 
standard model of cosmology, several "semi-public" or private codes, such as NUC123 [2] 
or PArthENoPE [3], already exist. 

We provide here a public code, AlterBBN, which enables the calculation of the abun- 
dances of the elements in the standard cosmological model, as well as in alternative cosmo- 
logical models such as quintessence or reheating scenarios. Indeed, many different phenom- 
ena in the Early Universe could have modified the Universe properties at the time of BBN, 
and in such cases the calculations could lead to abundances of the elements different from 
the predictions in the standard cosmological model. One of the main features of AlterBBN 
is that it provides calculations in such cosmological models, which can enable the user to 
test alternative scenarios and to constrain their parameters. 

The code AlterBBN has been inspired by NUC123, and we use similar calculation tech- 
niques throughout the program. In the following, we describe first the physics of Big-Bang 
nucleosynthesis, and present the set of equations in the standard model of cosmology. We 
then discuss alternative cosmological models, with a modified expansion rate or entropy 
production. We describe also the content of the AlterBBN package, as well as the list of 
its main routines. We will then explain how to compile and use AlterBBN, and the inputs 
and outputs of the program will be introduced. Finally, some examples of results obtained 
with AlterBBN will be given. In the Appendices, we give the list of the nuclear reactions 
used in AlterBBN, as well as experimental limits which can be considered to constrain the 
calculated abundances of the elements. 



2 BBN physics 

In the following section, we use the natural unit system c = h = k = 1. 



2.1 BBN in the Cosmological Standard Model 

All the nuclear reactions implemented in AlterBBN have evolution equations of the form: 

m M Z t + rij A iZj <— ► n k Ak Z k + n { Al Z { . (1) 
Therefore, the abundance change for nuclide i is given by [2]: 




where the nuclide abundance is Y{ = Xi/Ai, Xi being the mass fraction in nuclide % and 

ij 



Ai its atomic number. Uj is the number of nuclides i involved in the reaction, T^- is the 
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forward reaction rate and Tj k the reverse rate. 

The expansion rate of the Universe is given by the Friedmann equation: 



a\ 8vrG 



-Ptot 



a J 3 

and the equation of energy conservation can be written as: 

d 3 d 
Jt (ptota) + Ptot Jt 

where 
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We have 



Ptot = P-y + Pu + Pb+ (Pe- + Pe+) , 



Ptot = Pj + {P e - + Pe+) + Pb • 
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and in absence of neutrino degeneracy: 
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where N v = 3 is the number of neutrino families. The sum of the electron and positron 
densities is 
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and the sum of their pressures is 
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where z = m e /T, 



M(z) 



^ 4 (z) + ^ 2 (*) 



the being the modified Bessel functions. The charge conservation gives 

n e - - n e + = M A hT 3 S , 
where Ma is the Avogadro number, 

S = ZiYi , 

i 

and the difference between the electron and positron densities is 

2 °° 

= ^miy*(-l) n+l s\vh{ncb e )L{nz) , 
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(11) 

(12) 
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from which <j) e can be determined using: 

d(f) e d(f) e dT d<p e da d(j) e dS 

The baryon energy density reads 



Pb = hT 3 



1 + E 



( AM; 



+ {T)Yi 



where AMj is the mass excess of nuclide i, and the pressure is 



P b = hT 3 (^-CTY^Y^ , 



with C = 1-388 x 1(T 4 . h ~ p 6 /T 3 ~ l/a 3 /T 3 can be determined by the equation: 

dh 



dt 



-3h 



Ida 1 dT 

adi + f~dt 



(15) 



(16) 



(17) 



(18) 



To solve the system of differential equations which gives the abundances of the elements, 
we assume the following initial conditions: 



h( Ti ) = M U ^ V , 



(19) 



where M u is the unit atomic mass, n 7 is the number density of photons and r\ the baryon- 
to-photon ratio. Also, 



MTl) ~ ~2~ zf E,T=i(-i)" +1 ^(^) ' 

and the initial abundances of protons and neutrons are 



(20) 



(21) 



where q = m n — m p . The other quantities needed at this point are the initial baryon energy 
density 

Pb{T>i) « h(T t )Tf , (22) 

and the time 

U = (yi-KGof^T- 2 , (23) 
where a is the Stefan-Boltzmann constant. 



To solve this system of equations, a linearization of the nuclear abundance differential 
equations is performed, followed by a Runge-Kutta integration. 



2.2 Modified Cosmological Scenarios 

The Early Universe is a relatively unknown period, and therefore the properties of the 
Universe at that period are not well-known either. We consider in the following several 
alternative descriptions of the Early Universe properties. 
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2.2.1 Modified expansion rate 

Many phenomena at the time of BBN could have effects on the expansion rate of the 
Universe. In such cases, the Friedmann equation is affected, and we parametrize the modi- 
fication of the expansion by adding an effective dark energy density po to the total energy 
of the Universe. The Friedmann equation then reads: 

H 2 = I -J = —{p tot + p D ). (24) 

We adopt the parametrizations described in [HE] for pr,: 

pD = K p pt o t(T )l—J , (25) 

where Tq = 1 MeV. k p is therefore the ratio of effective dark energy density over total 
energy density, and n p is a parameter describing the behavior of this density. We refer the 
reader to |5 [5] for detailed descriptions and discussions of this parametrization. 



2.2.2 Modified entropy content 

The entropy content can also receive various contributions in the BBN period. We parametrize 
the contributions by considering an effective dark entropy density sd- The energy conser- 
vation equation reads in this case: 



- L {Ptota°) + Ptot-fXa ) - a° 



-T-(s D a 3 ) = 0. (26) 

T=cst dt 



dV r " u " ' w "dV ' dt 

We adopt two different parametrization for sd- The first one has been introduced in [5]: 

/ T ys 

s D = n s s rad (T ) 1 — 1 , (27) 



where Tq = 1 MeV and 

s rad {T) = h^(Ty-^T\ (28) 

^eff being the effective number of entropic degrees of freedom of radiation, k s the ratio of 
effective dark entropy density over radiation entropy density, and n s a parameter describ- 
ing the behavior of this density. We refer the reader to [5] for detailed descriptions and 
discussions of this parametrization. 

The second parametrization is inspired by reheating scenarios in which the entropy pro- 
duction T,£) evolves like [H] 

E D (T) = « s S rad (T ) f — J (29) 

where Tq = 1 MeV. is the ratio of effective dark entropy production over radiation en- 
tropy production, and ns is a parameter describing the behavior of this entropy production 
(ns ~ — 1 in most reheating scenarios). The radiation entropy production reads: 

/4vr 3 G \ 1/2 
S rad (T ) = ( — g eS (T )\ TQ 2 s rad (T ) , (30) 
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g e Q being the effective number of degrees of freedom of radiation, and where the entropy 
production stops at T = T r . The function £d(T) is related to the entropy density by: 

T 1/2, 



where 



- Wj^ 3 / *r , (si) 

Prad eff 



1/2 _ hefL ( Tdh^- 



This parametrization is further described in [6]. 



3 Content of the AlterBBN package 

AlterBBN is based on the library libbbn.a, which contains all the procedure necessary to 
the abundance computation, and is provided with five different programs. 



3.1 Parameter structure 

The main structure of the package AlterBBN is defined in src/ include. h and is provided 
below: 

typedef struct relicparam 
{ 

int entropy_model ; 

double ddO , ndd , Tdend ; 

double sdO,nsd,Tsend; 

double SigmadO , nSigmad , TSigmaend ; 

double ntO,nnt,Tnend; 

double mgravitino; 

double table.ef f [276] [3] ; 

double etaO; 

double nbnu; 

double lif e_neutron; 

> 

relicparam; 

This structure contains all the parameters which are necessary to the computation of the 
abundances of the elements in the standard cosmological model as well as in alternative 
models. 



3.2 Main routines 

The main routines of the libbbn.a library, which are needed for the abundance calculation, 
are detailed in the following: 
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• void Init_cosmomodel (struct relicparam* paramrelic) 



This function initializes the paramrelic structure, setting all the parameters to 0, 
apart from the number of neutrinos which is set to 3, and the baryon-to-photon ratio 
?? = 6.19 x 10" 10 pp. 



• void Init_cosmomodel_param(double eta, double nbnu, struct relicparam* paramrelic) 



This routine specifies values of the paramrelic structure, enabling the user to modify 
the baryon-to-photon ratio eta and the number of neutrinos nbnu. 



• void Init_dark_density (double ddO, double ndd, double T_end, struct relicparam* 
paramrelic) 



This routine changes values in the paramrelic structure, initializing the effective 
dark density parameters described in Appendix 2.2.1, with K p =ddO, n p =ndd, and 
where T_end is a temperature cut at which the dark density is set to 0. 



• void Init_dark_entropy (double sdO, double nsd, double T_end, struct relicparam* 
paramrelic) 



This routine changes values in the paramrelic structure, initializing the effective 
dark entropy density parameters described in Appendix 2.2.2, with K s =sdO, n s =nsd, 
and where T_end is a temperature cut at which the dark entropy density is set to 0. 



• void Init_dark_entropySigmaD (double SigmadO, double nSigmad, double T_end, 
struct relicparam* paramrelic) 



This routine specifies values of the paramrelic structure, initializing the effective 
dark entropy production parameters described in Appendix 2.2.2, with Ks=SigmadO, 
ns=nSigmad, and where T_end is a temperature cut at which the dark entropy pro- 
duction is set to 0. 



• void rate_weak (double f [] ) 



This procedure computes the decay rates of the weak interaction reactions, and stores 
them in f . 



• void rate_pn(struct relicparam* paramrelic, double f [] , double r[], double T9) 



This procedure computes the decay rate and reverse rate ofp ^ n reactions, and 
stores them in f and r respectively, at a temperature T9 (in unit of 10 9 K). 
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• void rate_all (double f [] , double T9) 

This procedure computes all the other reactions (summarized in Appendix [A]) at 
temperature T9 and stores them in f . 

• int nucl (struct relicparam* paramrelic, double *eta, double *H2_H, double *He3_H2, 
double *Yp, double *Li7_H, double *Li6_Li7) 

This function is the main routine of the program, as it calculates using the paramrelic 
structure the baryon-to-photon ratio eta, the ratio of the element abundances H2_H, 
He3_H2, Li7_H, Li6_Li7, and the helium abundance Y_p. It return if the calculation 
failed, or 1 otherwise. 

• int bbn_excluded (struct relicparam* paramrelic) 

This function is a "container" function which calls the routine nucl and returns 
if the abundances satisfy the constraints of Appendix [Bj 1 otherwise, and -1 if the 
calculation failed. 

4 Compilation and installation instructions 

AlterBBN is written for a C compiler respecting the C99 standard, and it has been tested 
successfully with the GNU C and the Intel C Compilers on Linux and Mac. The package 
can be downloaded from: 



http : / / superiso . in2p3 . f r/relic/ alterbbn 



After unpacking, the following main directory is created: 
alterbbn_vX . X 

This directory contains the src/ directory, in which all the source files can be found. The 

main directory contains also a Makefile, a README, five sample main programs (stand_cosmo . c, 

alter_eta.c, alter_neutrinos . c, alter_standmod. c and alter_reheating. c). The 

compilation options should be defined in the Makefile, and in particular the C compiler 

name and its specific flags, if needed. 

Additional information can be found in the README file. 

To compile the library, type 

make 

This creates libbbn.a in src/. Then, to compile one of the five programs provided in 
the main directory, type 
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make name or make name . c 



where name can be stand_cosmo, alter_eta, alter_neutrinos, alter_standmod or alter_reheating. 
This generates an executable program with the .x extension: 

stand_cosmo .x calculates the abundance of the elements in the standard cosmological 
model, assuming the baryon-to-photon ratio measured by WMAP. 

alter_eta.x computes the abundance of the elements in the standard cosmological model, 
with the baryon-to-photon ratio given by the user. 

alter_neutrinos .x calculates the abundance of the elements in the standard cosmological 
model, and give the user the possibility to alter the number of neutrinos. 

alter_standmod.x and alter_reheating computes the abundance of the elements in cos- 
mological scenarios with modified expansion rates and entropy contents. 

5 Input and output description 

In the following, we describe the input and output of the main programs. 

5.1 Standard cosmology 

The program stand_cosmo .x computes the abundance of the elements in the standard 
cosmological model assuming the baryon-to-photon ratio r] = 6.19 x 10~ 10 measured by 
WMAP [I]. No parameter is needed by this program. Running the program with: 
. /stand_cosmo . x 
returns 

Yp H2/H He3/H2 Li7/H Li6/Li7 

2.448e-01 2.544e-05 4.079e-01 4.477e-10 2.429e-05 

Compatible with BBN constraints 

5.2 Standard cosmology with r] modification 

The program alter_eta.x computes the abundance of the elements in the standard cos- 
mological model, and needs the baryon-to-photon ratio r\ as input. Running the program 
with: 

. /alter_eta.x 3e-10 

returns 

Yp H2/H He3/H2 Li7/H Li6/Li7 

2.374e-01 7.967e-05 2.020e-01 1.168e-10 2.692e-04 

Excluded by BBN constraints 

5.3 Standard cosmology with modified neutrino number 

The program alter_neutrinos .x computes the abundance of the elements in the standard 
cosmological model with a modified neutrino number. It needs as input the number of 
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neutrino families N v . 
Running the program with: 
. /alter_neutrinos .x 4.5 

returns 

Yp H2/H He3/H2 Li7/H Li6/Li7 

2.629e-01 3.053e-05 3.614e-01 3.914e-10 3.641e-05 

Excluded by BBN constraints 



5.4 Effective energy and entropy densities 

The program alter_standmod.x computes the abundance of the elements while adding to 
the standard cosmological model an effective energy density such that 

PD = K p p T ad{TBBN)(T /T B BN) n " , (33) 

and/or an effective entropy density 

s D = K s s rad (T B BN)(T/T B BN) ns , (34) 

which modify the Early Universe properties without having observational consequences if 
chosen adequately [H [5]. A description of the model and of the related equations can be 
found in Sections 2.2.1 and 2.2.2 The necessary arguments to this program are: 

• K p : ratio of dark energy density over radiation energy density at BBN time, 

• n p : dark energy density decrease exponent, 

• k s : ratio of dark entropy density over radiation entropy density at BBN time, 

• n s : dark entropy density decrease exponent. 

Optional arguments can also be given: 

• T p : temperature in GeV below which the dark energy density is set to 0, 

• T s : temperature in GeV below which the dark entropy density is set to 0. 

For example, the command . /alter_standmod.x 0.1 6 l.e-3 5 

returns 

Yp H2/H He3/H2 Li7/H Li6/Li7 

2.315e-01 4.558e-05 2.831e-01 1.933e-10 9.173e-05 

Compatible with BBN constraints 



5.5 Effective reheating model 

The program alter_standmod. x computes the abundance of the elements while adding to 
the standard cosmological model an effective energy density such that 

PD = K p p rad (T B BN){T '/Tbbn)"" , (35) 
and/or an effective entropy production 

= HT^rad{TBBN){T /Tbbn)^ , (36) 
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which modify the Early Universe properties without having observational consequences if 
chosen adequately [6 J • A description of the model and of the related equations can be found 



in Sections 2.2.1 and 2.2.2 Note that is related to sd by Eq. (31). The necessary 



arguments to this program are: 

• K p : ratio of dark energy density over radiation energy density at BBN time, 

• n p : dark energy density decrease exponent, 

• ky,'- ratio of dark entropy production over radiation entropy production at BBN time, 

• n-£,: dark entropy production exponent, 

• T r : temperature in GeV below which the dark energy density and the entropy pro- 
duction are set to 0. 

For example, the command ./alter_reheating.x 0.1 6 0.1 -1 le-3 

returns 

Yp H2/H He3/H2 Li7/H Li6/Li7 

2.561e-01 2.615e-05 4.007e-01 4.590e-10 2.583e-05 

Compatible with BBN constraints 

In order to perform scans in different cosmological models, the user is invited to write 
his/her own programs using the above main programs as guidelines. 



6 Results 

AlterBBN has been thoroughly tested, and provides results in agreement with those of 
NUC123 H]. 

Considering the different cosmological models implemented, it is possible to perform 
scans on the cosmological parameters of these models to determine constraints. For example, 
in Fig. [TJ the current limits on the modified expansion properties from the Y p and 2 H/H 
BBN constraints are presented. The area between the back lines on the left plot and the 
area on the top of the black lines on the right plot lead to unfavored element abundances. 
The constraints of Appendix [B] are used in this figure. 

Also, AlterBBN will soon be interfaced and provided in the Superlso Relic package [7], 
so that the implemented cosmological models can be tested at the same time by BBN 
constraints as well as by particle physics constraints. 
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Figure 1: Constraints from Y p (left) and 2 H/H (right) on the effective dark energy density 
parameters (n p , n p ). The parameter regions excluded by BBN are located between the black 
lines for Y p , and in the upper left corner for 2 H/H. The colours correspond to different 
values of Y p and 2 H/H. 
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Appendix A Nuclear reaction network 

AlterBBN currently includes a network of 88 nuclear reactions, which are gathered in the 
following table. 
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Table 1: Network of nuclear reactions implemented in AlterBBN. 
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Appendix B BBN constraints 

The following conservative constraints |16| are used in the function bbn_excluded: 

0.240 < Y p < 0.258 , 1.2 x 10~ 5 < 2 H/H < 5.3 x 10~ 5 , (37) 
0.57 < 3 H/ 2 H < 1.52 , 7 Li/H > 0.85 x 10^ 10 , 6 Li/ 7 Li < 0.66 , 

for the helium abundance Y p and the primordial 2 H/H, 3 H/ 2 H, 'Li/H and e Li/ 7 Li ratios. 
It is possible for the user to change these constraints in routine bbn_excluded which can 
be found in src/bbn.c. 
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